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Foil  angle  relative  to  pod 
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ABSTRACT 


As  part  of  an  investigation  into  minimizing  control 
power  requirements  for  hydrofoil  craft,  a technique  has 
been  developed  for  obtaining  complete  hydrodynamic  lift 
and  moment  characteristics  of  a flapped  hydrofoil  from 
equilibrium  data  obtained  in  flap  incidence  control  experi- 
ments with  a freely  pivoting  foil.  This  method  employs  a 
simple  mathematical  model  and  various  crossplots  of  the 
data  to  separate  the  effects  of  camber,  foil  incidence 
angle,  and  flap  angle  using  a minimum  amount  of  data. 

Existing  model  data  for  the  AGEH-1  main  strut -pod -foil 
system  under  subcavitating  conditions  have  been  analyzed 
by  this  method,  and  effects  of  depth  and  velocity  are 

described.  Velocity  effects  were  attributed  to  air  and/or 

} ► 

water  flow  out  of  the  pod.  Relatively  high  values  of  flap 
lift  effectiveness  (0.3  to  0.4)  were  obtained.  Centers  of 
pressure  (lines  of  action)  for  lift  due  to  camber,  angle 
of  attack,  and  flap  angle  are  presented.  Recommendations 
are  given  for  improving  the  accuracy  of  the  measurement 
and  analysis  techniques. 


ADMINISTRATIVE  INFORMATION 


m 


The  work  presented  here  was  sponsored  by  the  Naval  Sea  Systems 
Command  (NAVSEA  SYSCOM)  under  Project  S4606,  Task  01703,  and  performed 
under  the  David  W.  Taylor  Naval  Ship  Research  and  Development  Center 
(DTNSRDC)  Work  Unit  1153-503. 

INTRODUCTION 

Of  recent  interest  to  the  U.S.  Navy  have  been  methods  of  reducing 
the  control  power  requirements  for  hydrofoil  craft.  In  a theoretical 
study  of  four  different  hydrofoil  configurations  in  a seaway.  Remington 
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Hi 
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and  Bender*  concluded  that  the  configuration  using  a trailing  edge 

flap  for  incidence  control  of  a freely  pivoting  foil  required  the 

least  amount  of  power.  An  existing  model  of  the  AGEH  main  strut -pod- 

foil  with  a trailing  edge  flap  on  the  foil  has  been  modified  for 

flap  incidence  control . The  foil  was  made  freely  pivoting  in  pitch 

(incidence  angle)  about  the  existing  transverse  axis  through  the  pod. 

The  motion  of  the  flap  and  foil  angles  were  coupled  together  by  a 

series  of  linkages  to  obtain  a passive  control  system.  An  experiment 

performed  at  DTNSRDC  using  this  modified  model  has  been  reported  by 
2 

Coder  and  Wisler.  The  purpose  of  that  experiment,  labeled  here  as 
the  "1975  Flap  Incidence  Control  Experiment,"  was  to  determine  whether 
equilibrium  flap-foil  angle  conditions  ("stable  flying  situations") 
could  be  obtained.  The  velocities  were  low  (7.31  to  12.86  m/s)  so  that 
the  foil  was  fully  wetted  (no  cavitation)  for  most  of  the  conditions 
tested.  Equilibrium  conditions  were  obtained  and  are  reported  in 
Reference  2.  While  those  data  were  being  reduced,  it  was  discerned 
that  the  moment  at  zero  flap  and  foil  angles  and  the  variance  of  moment 
with  flap  angle  alone  could  be  obtained  from  the  data  using  what  will 
be  called  here  the  "indirect  method."  After  those  results  were  published, 
this  "indirect  method"  was  expanded  to  include  consideration  of  the  lift 
data  (not  published  in  Reference  2 but  included  here  in  Appendix  A)  to 
determine  other  hydrodynamic  loading  characteristics. 


References  are  listed  on  page  79. 

This  model  was  previously  used  for  work  reported  informally  by 
A.C.  Conolly  in  DTNSRDC  Report  SPD-332-H-02  in  February  1973. 
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This  "indirect  method"  is  based  on  that  two-dimensional  repre- 
sentation of  the  foil  developed  in  Coder,  et  al."*  - the  lift  and 
moment  were  assumed  to  be  the  sum  total  of  the  parts  due  to  foil  camber, 
foil  angle  of  attack,  and  tab/flap.  Here  the  equations  for  lift  and 
moment  are  written  in  coefficient  form  as  given  in  Equations  [1]  and  [2] 

of  the  ANALYSIS  section.  Using  the  straight-line  fits  to  the  flap 

2 

angle-foil  angle  equilibrium  data  from  Coder  and  Wisler,  and  appro- 
priate crossplots  of  those  data  and  the  lift  data  from  Appendix  A,  the 
six  coefficients  on  the  right-hand  side  of  Equations  [1]  and  [2]  are 
determined.  Also,  it  will  be  seen  that  the  lift  and  moment  effective- 
ness and  the  lines  of  action  of  the  various  lift  components  may  be 
determined  from  the  resulting  coefficients. 

The  objective  of  the  present  work  is  to  determine  with  actual  data 
whether  this  "indirect  method"  produces  reasonable  results.  Since  it 
is  planned  to  conduct  further  flap  incidence  control  experiments  with  an 
improved  model  and  under  cavitating  conditions,  this  "indirect  method" 
may  prove  to  be  a viable  method  to  determine  these  lift  and  moment 
coefficients  as  a byproduct  of  studying  the  equilibrium  conditions  and 
the  control  system.  Further,  it  may  prove  to  be  more  advantageous  from 
a testing  point  of  view  to  determine  the  loading  characteristics  using 
this  "indirect  method"  since  much  data  is  obtained  rather  quickly  and 
efficiently. 

The  remainder  of  the  report  describes  this  "indirect  method," 
presents  and  discusses  the  results  for  the  data  from  the  1975  Flap 
Incidence  Control  Experiments,  and  makes  recommendations  as  to  the  future 
employment  of  this  method. 

ANALYSIS 

The  lift  and  moment  in  coefficient  form  for  a fully  wetted  flapped 
foil  at  small  angles  may  be  represented  by 
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where 


C 


L 


q 

s 

U 

oc 


[1] 

St/  * CHC  * CH?f 

[2] 

is  lift  coefficient  («  lift/qS) 

is  foil  lift-curve  slope  (=30^/3 o£) 

is  camber  lift  coefficient  (also  includes  thickness  effects) 

is  flap  lift-curve  slope  (=  3Cl/3<9) 

is  moment  coefficient  (=  moment/q  Sc) 

is  foil  moment-curve  slope  (=3C„/3«) 

is  camber  moment  coefficient 

is  flap  moment-curve  slope  (=3CH/3^3) 

is  mean  geometric  chord  based  on  foil  planform  extended 
to  pod  centerline 

2 

is  dynamic  pressure  (=  (1/2)^  U ) 
is  planform  area  of  foil  extended  to  pod  centerline 
is  free  stream  velocity 
is  foil  angle 


^ is  flap  angle 

^ is  density  of  fluid 

The  data  reported  by  Coder  and  Wisler2  pertain  to  a flapped,  freely- 
pivoting  foil  and  are  presented  in  the  form  of  foil  angle  versus  flap 
angle  to  maintain  hydrodynamic -external  bias  moment  equilibrium,  i.e., 

4 


in  the  form 


0 = m OC  + b [3] 

where  m is  the  slope  (=  ?(»/2oc) 
b is  the  intercept  (= 

From  these  data  alone,  CH(,  and  CH^  were  calculated  and  are  presented 
in  that  report.  The  present  report  also  presents  (in  Appendix  A)  lift 
coefficient  (C^)  versus  foil  angle  (.OC)  for  that  experiment.  The 
slopes  and  intercepts  of  average  straight-lines  drawn  through  these  data 

are  presented  in  Table  1 together  with  pertinent  data  from  Coder  and 

2 

Wisler.  From  all  these  data,  the  coefficients  in  Equations  [1]  and 
[2]  and  centers  of  pressure  (lines  of  action)  for  the  various  lift  com- 
ponents may  be  determined.  A schematic  showing  the  analysis  technique  is 
shown  in  Figure  1 and  will  now  be  described. 

LIFT  COEFFICIENTS 

If  Equation  [3]  is  used  with  Equation  [1] 


Thus  we  have 


CL  = CLof°d  + CLC  + CL  (3  * + b') 

" (CL*  * (CLC  + b(V 


(ClV  = 0 = CLC  + CL£  b 


dCL/d*  = CL«  + CL£  m 


[4] 


[5] 


In  Figure  2,  (C^)^  _ Q from  Table  1 is  plotted  versus  b from  Table  1 from 
which  C^  and  C^  were  determined  as  the  slope  and  intercept,  respectively 
and  listed  in  Table  2.  In  Figure  3,  6CL/du  from  Table  1 is  plotted 
versus  m from  Table  1.  As  seen  in  the  figure,  it  is  impossible  to 
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LIFT  COEFFICIENT  FOR  ZERO  FOIL  ANGLE,  (C,  L,  n (FROM  TABLE  1) 


FLAP  ANGLE  IN  DEG  FOR  ZERO  FOIL  ANGLE,  (01,  « OR  b (FROM  TABLE  1) 

«*•  u 


FIGURE  2 - GRAPHICAL  DETERMINATION  OF  LIFT  AT  ZERO  FLAP 
AND  FOIL  ANGLES  (CLC)  AND  FLAP-LIFT-CURVE  SLOPE  (C^) 
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COEFFICIENT  FOR  ZERO  FOIL  ANGLE,  (C,  V*  A (FROM  TABLE  1) 


WKmmm 


h/c  * 0.578 


slope  * cLg  * (^/se)^  0 (determined  here) 


INTERCEPT  = C, 


(ClWb=  0 (DETERMINED  HERE) 


FLAP  ANGLE  IN  DEG  FOR  ZERO  FOIL  ANGLE,  (6)^  Q OR  b (FROM  TABLE  1) 
FIGURE  2b  - VELOCITY  « 9.00  m/s  (17.5  KNOTS) 


LIFT  COEFFICIENT  FOR  ZERO  FOIL  ANGLE,  (C^^  (FROM  TABLE  1) 
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FIGURE  3a  - VELOCITY  =7.31  m/s  (14.2  KNOTS) 


38/Mt  OR  m IN  DEG/DEG  (FROM  TABLE  1) 
FIGURE  3b  - VELOCITY  = 9.00  m/s  (17.5  KNOTS) 

FIGURE  3 - GRAPHICAL  DETERMINATION  OF  FOIL-LIFT-CURVE  SLOPE  (Cj^) 
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3C  /3«  IN  DEG'1  (FROM  TABLE  1)  3C  /3*  IN  DEG'  (FROM  TABLE  1) 


draw  a straight  line  through  the  data  so  that  an  accurate  value  for 
and  CLo(  could  be  obtained  from  the  slope  and  intercept,  respectively. 
However,  since  the  appropriate  slope  for  the  straight  line,  C.  _ , has 

L p 

already  been  determined  from  Figure  2,  it  was  used  to  draw  the  straight 
line  in  Figure  3 so  that  the  intercept,  C.  , could  be  obtained.  C.  . 

L|A  LX 

determined  this  way  is  also  listed  in  Table  2. 

MOMENT  COEFFICIENTS 

Since  the  hydrodynamic  moment  (C„)  and  the  bias  moment  (Cw)  together 

H M 

are  zero  for  a free  foil. 


CH  * -CM 


From  Equation  [2], 


'CM  = + CHC  + CH^  f 


$ = ’^CM  + CHC^/CHp  ^ 


From  Equation  [3]  we  have 


m CC  + b 


m = "(CH*/CH^  5 
b " -fCM  + CHC^CH( 


CM  = ^HC5  + (_CH^  } b 

- 

CH*  = "m  CH^ 

By  plotting  C„  from  Table  1 versus  b from  Table  1 in  Figure  4,  -C,,- 

M Hfl 

and  -CHC  were  determined  as  the  slope  and  intercept,  respectively.  Then 
using  m averaged  over  the  bias  moments  from  Table  1 and  C„^  as  determined 


BIAS  MOMENT  COEFFICIENT,  C (FROM  TABLE  1) 


SLOPE  » Cne  - 0 (DETERMINED  HERE) 

_j j ^ — -/H { \ - { / 

INTERCEPT  * C = (CJ-  0 n (DETERMINED  HERE) 


+ 12.86 


-6.0  -4.0 

FLAP  ANGLE  IN  DEG  FOR  ZERO  FOIL  AN 

FIGURE  4a  - h/~  ^78 


Xa  0 OR  b (FROM  TABLE  1) 


FIGURE  4 - GRAPHICAL  DETERMINATION  OF  MOMENT  AT  ZERO  FLAP 
AND  FOIL  ANGLES  (C^)  AND  FLAP -MOMENT- CURVE  SLOPE  (C^) 


BIAS  MOMENT  COEFFICIENT,  Cu  (FROM  TABLE  1) 


from  Figure  4,  was  determined.  These  values  for  moment  coefficients 

are  also  listed  in  Table  2. 


CENTERS  OF  PRESSURE  OF  LIFT  (LINES  OF  ACTION  FOR  THE  THREE  LIFT  COMPONENTS) 


Using  the  above  sketch,  the  moment  equation.  Equation  [2],  may  also 
be  written  as 


where 


X 


CH  ' CI**  (X-Y)  * CLC 


(x-z)  . CL^  (9  (X-W) 


is  the  foil  pivot  location  from  the  leading  edge  in  mean 
geometric  chords, 


Y is  the  planform  lift  center  of  pressure  location  from  the 
leading  edge  in  mean  geometric  chords, 

Z is  the  camber  lift  center  of  pressure  location  from  the 
leading  edge  in  mean  geometric  chords, 

W is  the  flap  lift  center  of  pressure  location  from  the 
leading  edge  in  mean  geometric  chords 


Thus, 


“H* 


“HC 


‘ CL*  <X"Y> 
CLC  (X'Z> 


CH£  * CL£  (X‘W) 
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RESULTS 

When  viewing  the  results  of  the  analysis  presented  in  Figures  5 
through  8,  it  should  be  noted  that  the  analyzed  results  for  12.86  m/s  are 
not  reliable  since  data  for  only  two  bias  moments  (instead  of  four) 
were  obtained  and  analyzed.  For  this  reason,  the  hand- faired  curves 
through  the  analyzed  results  in  the  figures  are  dashed  when  extending 
the  results  from  10.55  m/s  to  12.86  m/s.  For  most  cases,  these  results 
for  12.86  m/s,  do  however,  follow  the  general  trends  established  from 
the  lower  velocities. 

As  shown  in  Figure  5,  over  the  range  of  model  depths  and  velocities, 
the  flap  lift-curve  slope  (Cj^  ^ increases  slightly  with  increasing  depth 
and  decreases  slightly  with  increasing  velocity.  The  foil  lift -curve 
slope  (Cj^)  increases  with  increasing  depth  and  increasing  velocity. 
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LIFT  CURVE  SLOPE  IN  RAD  AND  LIFT  EFFEC1 


FIGURE  S - LIFT  CURVE  SLOPES  AND  FLAP  LIFT  EFFECTIVENESS  VERSUS  VELOCITY 


MOMENT  CURVE  SLOPE  IN  RAD  AND  MOMENT  EFFECTIVENE! 


VELOCITY  IN  METERS/SECOND 


FIGURE  6 - MOMENT  CURVE  SLOPES  AND  FLAP  MOMENT  EFFECTIVENESS  VERSUS  VELOCITY 


LIFT  AND  MOMENT  COEFFICIENTS 


■HMMIliiiMHWHMlH 


O—'Ny 6.0  8.0  10.0  12.0  14.0 

VELOCITY  IN  METERS/ SECOND 


FIGURE  7 - LIFT  AND  MOMENT  FOR  ZERO  FLAP  AND  FOIL  ANGLES  VERSUS  VELOCITY 
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0 — — 6.0  8.0  10.0  12.0  14.0 

VELOCITY  IN  METERS/SECOND 

FIGURE  8a  - DUE  TO  PLANFORM  (Y)  AND  CAMBER  (Z) 

FIGURE  8 - CENTER  OF  PRESSURE  OF  LIFT  DUE  TO 
PLANFORM  00.  CAMBER  (Z) , AND  FLAP  (W) 
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DISTANCE  TO  CENTER  OF  PRESSURE  IN  MEAN  GEOMETRIC  CHORDS 


The  resulting  flap  lift  effectiveness  (E^  = /C^)  shows  no  depth 

dependence  but  a slight  decrease  with  increasing  velocity. 

In  Figure  6,  depth  has  a very  small  effect  on  the  flap  moment- 
curve  slope  (CH^  ) and  the  foil  moment-curve  slope  Both  "Ch^ 

and  increase  significantly  with  increasing  velocity.  When  the 
flap  moment  effectiveness  (E^  = C^/C^)  *s  calcu^ate^»  the  effect 
of  depth  is  more  easily  discernible  --  its  magnitude  is  larger  for  larger 
depths.  The  magnitude  of  EH  decreases  slightly  with  increasing  velo- 
city. 

The  lift  and  moment  coefficients  for  zero  flap  and  foil  angles  are 
shown  in  Figure  7.  The  magnitudes  of  the  coefficients  are  larger  for  the 
larger  depth.  The  lift  coefficient  for  zero  flap  and  foil  angles  (C^) 
decreases  somewhat  and  the  negative  of  the  moment  coefficient  for  zero 
flap  and  foil  angle  (-CHC)  increases  significantly  with  increasing 
velocity.  This  large  increase  in  the  magnitude  of  this  moment  coefficient 
is  due  to  the  aft  movement  of  the  center  of  pressure  of  lift  due  to  camber 
(Z  * X as  shown  in  Figure  8.  This  center  of  pressure  is  slightly 

more  aft  for  the  larger  depth. 

The  center  of  pressure  of  lift  due  to  planform  (Y  = X - C„  . /C.^  ) 

HCC  LAC 

moves  slightly  forward  with  velocity  up  to  10.55  m/s  and  approaches  the 
leading  edge  (near  zero)  for  12.86  m/s.  The  center  of  pressure  of  lift 
due  to  flap  (W  = X -C^  5 slightly  more  aft  at  the  larger  depth 

and  moves  aft  with  increasing  velocity.  Thus,  it  is  seen  that  the  two 
centers  of  pressure  Z and  W shift  with  velocity  such  as  to  cause  a larger 
negative  moment  (pitching  leading  edge  downward)  with  increasing  velocity 
and  Y shifts  slightly  more  forward  to  cause  a larger  positive  moment. 


DISCUSSION 

The  foregoing  analysis  and  results  show  that  it  is  possible  to 
obtain  the  six  hydrodynamic  lift  and  moment  coefficients  from  a foil 
that  is  freely  pivoting,  has  a controllable  flap  with  suitable  feedback 
control  system,  and  has  provision  for  externally  applied  bias  moments. 
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The  coefficients  such  obtained  are  not  as  accurate  as  those  that  might 
be  obtained  by  direct  measurement.  For  example,  the  flap  lift  curve  slope 
could  be  obtained  more  accurately  by  fixing  the  foil  at  a set  angle  (say 
zero)  and  measuring  the  lift  while  varying  the  flap  angle.  The  present 
method  does,  however,  provide  a simple  determination  of  the  coefficients 
from  a relatively  small  body  of  data.  When  using  this  method,  it  would 
be  best  to  ensure  that  for  every  condition  (velocity  and  depth)  data  were 
obtained  for  at  least  three  bias  moments  and  preferably  four. 

The  most  significant  result  of  the  analyzed  results  is  the  drastic 

shift  in  centers  of  pressure  due  to  velocity  (shown  in  Figure  8).  For 

a similar  foil  (but  without  a flap)  at  the  lower  velocities  (7.31  and 

3 

10.55  m/s),  it  was  determined  by  Coder,  et  al  that  Y (Planform  Center) 
was  near  the  quarter  chord  and  Z (Camber  Center)  was  at  about  the  three- 
quarter  chord  with  no  apparent  effect  of  velocity.  These  results  are 
shown  on  Figure  8a  for  a comparison  with  the  present  data.  The  present 
data  shows  a greater  dependence  on  velocity.  It  is  thought  that  the 
effect  of  velocity  on  the  centers  of  pressure  of  the  data  analyzed  here 
might  well  be  due  to  small  amounts  of  water  and/or  air  leaking  out  of  the 
pod  into  the  flow.  The  experiment  described  in  Reference  4 had  one  push- 
rod  passing  vertically  through  the  strut  which  was  a possible  path  for 
atmospheric  air  leakage  to  the  pod  and  foil.  An  attempt  was  made  to 
seal  the  strut  to  prevent  this  source  of  leakage.  Whenever  this  was 
done,  it  was  found  that  the  resulting  hydrodynamic  coefficients  became 
constant  for  the  two  velocities.  In  the  present  experiment,  discussed 
in  Reference  2,  it  was  more  difficult  to  seal  the  strut  because  two 
linkages  passed  vertically  through  the  strut  and  it  was  imperative  that 
the  sealing  technique  not  increase  friction  in  the  linkages.  Little,  if 
any,  air  was  observed  during  testing.  However,  it  is  apparent  now  that 
there  might  have  been  enough  leakage  to  influence  the  results. 

The  large  shifts  in  centers  of  pressure  at  the  higher  velocities 
(10.55  and  12.86  m/s)  may  be  partly  due  to  the  license  used  in  drawing 
the  straight  line  fits  through  the  data  during  the  analysis  procedure. 

For  example,  the  straight  lines  drawn  on  Figure  4 were  drawn  through  the 
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average  of  0 for  CM  ■ 0 of  the  three  lower  velocities.  It  is  seen 

that  the  resulting  lines  drawn  for  the  larger  velocities  of  10.55  and 
12.86  m/s  are  steeper  and  have  larger  intercepts  than  those  that  might  be 
drawn  through  the  actual  data  points  for  the  particular  velocity  alone. 

This  would  result  in  a larger  -C^,  -C^  and  CHflJ  (since  this  is  determined 
as  -m  than  what  ordinarily  might  be  determined.  Thus,  Y,  Z,  W,  and 

are  affected  by  this  interpretation  --  Y might  be  too  small,  and  Z,  W,  and 
-Eh  too  large.  This  would  tend  to  make  the  results  show  a weaker  dependence 
on  velocity. 

A detailed  error  analysis  is  performed  in  Appendix  B.  Estimated 
maximum  errors  in  the  fundamental  experimental  measurements  were  determined 
and  followed  through  the  analysis  step  by  step.  Error  bands  for  the 
hydrodynamic  coefficients,  centers  of  pressure,  and  effectiveness  were  thus 
obtained.  This  error  analysis  shows  that  there  could  be  considerable  errors 
apparently  resulting  mainly  from  the  inaccuracy  of  the  flap  angle  measurement. 
Improvement  in  this  measurement  and  perhaps  foil  angle  and  bias  moment  could 
improve  the  results  considerably. 

CONCLUSIONS 

1.  The  technique  developed  in  this  report  permits  determination  of 
the  hydrodynamic  loading  coefficients  characteristics  on  a freely  pivoting 
foil  with  controllable  flap,  a suitable  control  system,  and  external 
applied  bias  moments.  Measurement  of  lift,  foil  angle,  and  flap  angle  must 
be  made  and  the  bias  moment  known. 

2.  To  use  this  technique  effectively,  data  for  several  values  of  bias 
moment  are  required.  At  least  three  values,  and  preferably  four,  for  each 
velocity  and  depth  should  be  used. 

3.  A critical  step  in  the  analysis  is  drawing  the  straight  line  fits. 
Varying  the  approach  can  alter  the  results. 

4.  Error  bands  for  the  hydrodynamic  coefficient,  centers  of  pressure, 
and  effectivenesses  could  be  considerably  improved  by  an  improvement  in  the 
flap  angle  measurement.  Lesser  amounts  of  improvement  could  be  obtained  by 
improving  the  foil  angle  and  bias  moment  measurements. 
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5.  The  data  of  Reference  2 is  not  a good  representation  of  the 
full-scale  AGEH  system.  The  present  analysis  suggests  that  there  was  air/ 
water  leaking  down  the  strut  and  out  of  the  pod  into  the  flow. 

RECOMMENDATIONS 

1.  In  future  experiments  employing  the  present  technique,  data  for 
four  values  of  bias  moment  should  be  obtained  for  each  velocity  and  depth. 

2.  In  future  experiments  to  determine  hydrodynamic  coefficients,  etc., 
by  the  "indirect  method",  the  flap  angle  measurement  should  definitely  be 
improved.  Also  further  improvement  of  the  results  could  be  obtained  by 
improving  the  foil  angle  and  bias  moment  measurements. 

3.  The  experiment  should  be  extended  into  the  region  of  cavitating 
condition  to  see  what  effect  cavitation  has  on  the  hydrodynamic  character- 
istics and  the  "indirect  method"  analysis  technique. 
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APPENDIX  A 


LIFT  COEFFICIENT  VERSUS  FOIL  ANGLE  FROM  THE 
1975  FLAP  INCIDENCE  CONTROL  EXPERIMENT 
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APPENDIX  B 

ERROR  ANALYSIS  OF  "INDIRECT  METHOD" 

The  "indirect  method'bf  determining  the  hydrodynamic  coefficients, 
lift-  and  moment-curve  slopes,  centers  of  pressure,  and  flap  effective- 
ness is  shown  in  Figure  B1  (also  shown  as  Figure  1 of  main  report).  This 
analysis  is  based  on  the  linear  equations  developed  in  the  report  (Equa- 
tions [1]  through  [3]).  No  attempt  is  made  here  to  determine  the  errors 
introduced  by  the  basic  assumptions  of  those  equations,  but  rather  to 
follow  the  errors  due  to  measurement  of  the  fundamental  experimental 
parameters  through  the  analysis  procedure. 

The  estimated  maximum  errors *in  the  fundamental  experimental 
measurements  are  listed  in  Table  Bl.  Based  on  these  errors,  the  estimated 
maximum  errors  in  the  derived  parameter  may  be  determined  as  follows: 

1.  LOCATION  OF  PIVOT  POINT  (X) 

X = (Xc)/c 

dx/x  = d(xc)/(xc)  - dc/c” 

jdx/xj  4.  jd(xc)/(xc)j  + jdc/c| 

From  Table  Bl, 

|dx/x|  0.00036  + 0.00012  = 0.00048 
Thus,  maximum  percent  error  in  X is  about  0.05  percent. 


* These  estimated  maximum  errors  include  both  accuracy  and  precision 
errors. 
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2.  BIAS  MOMENT  COEFFICIENT  (CM) 

CM  = M/(1/2^U2Sc) 

__ 

dCM/CM  = dM/M  - df/f  -2dU/U  -dS/S  -dc/c 

|dCM/cM  | ^ |dM/M|  + |d^|  + 2 |dU/u|  ♦ |dS/s|  ♦ |dc/c  | 

From  Table  Bl, 


dS/CM  I - 


+ 0.00008  + 2 


Not  Defined  for  M = 

0 ft-lb 

0.14 

for  M = 

10  ft-lb 

0.07 

for  M = 

20  ft-lb 

o.os 

for  M = 

30  ft-lb 

0.0007 

for 

U = 14.2 

Kts 

0.0006 

for 

U = 17.5 

Kts 

0.0005 

for 

U = 20.5 

Kts 

0.0004 

for 

U = 25  Kts 

+ 0.00020  + 0.00012 


All  other  terms  on  the  right  side  of  the  equation  are  negligible  compared 
to  the  first  term,  so 

|dCM  | f |dM/M | x j CM  | = |dM/(l/2^U2  Sc)| 

The  maximum  error  in  Cw  is  calculated*  as  a function  of  M and  U and 

M 

is  shown  in  Table  B2. 


Using  1/2(0  Sc  » 1.090  lb-sec  /ft 
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TABLE  B2 

ESTIMATED  MAXIMUM  ERROR  IN  BIAS  MOMENT  COEFFICIENT  (CJ 


(ft-lb) 

0 

10 

20 

30 

14.2 

0.0022 

0.0023 

0.0022 

0.0025 

17.5 

0.0015 

0.0015 

0.0015 

O.OO06 

20.5 

0.0011 

0.0011 

0.0011 

0.0012 

25.0 

0.0007 

0.0007 

0.0007 

0.0008 

The  range  for  lift  is  0 ^ L 400  lbs,  so  the  first  term  on  the  right  is 
0.005  ^ 2/  L ^ oC 


All  other  terms  are  small  compared  to  this  one  so  we  may  write 
| dCL/CL | * ldL/Ll 

|dCL  | A JdL/L  | x |CL  I - |dL/(l/2  £ U2S)| 


The  maximum  error  in  is  calculated*  for  the  various  velocities  as 


* Using  1/2  ^ S = 1.512  lb-sec2/ft2 
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|dCL|  = 2/(1.512  U2(ft/sec)2) 

{0.0021  for  U = 14.2  Kts 
0.0015  for  U = 17.0  Kts 
0.0011  for  U = 20.5  Kts 
0.0007  for  U = 25.0  Kts 

4.  FOIL  ANGLE  (*) 

The  estimated  maximum  error  in  foil  angle  from  Table  B1  is  0.13  deg. 

5.  FLAP  ANGLE  (^) 

The  estimated  maximum  error  in  flap  angle  from  Table  B1  is  0.78  deg. 

From  the  derived  parameters  onward,  the  analysis  technique  involves 
various  calculations  and  plots  as  shown  in  Figure  Bl.  The  errors  in  the 
slopes  and  intercepts  obtained  from  the  plots  are  determined  as  followed. 

For  example,  take  the  general  case  of  a plot  of  one  variable  (F)  versus 
another  variable  (G)  as  shown  in  Figure  B2a.  A "best”  straight  line  is 
drawn  by  hand.  From  this  line,  "best"  values  for  slope  and  intercept  are 
determined.  These  are  the  values  listed  in  the  report.  Since  the 
variables  F and  G may  have  error,  the  data  points  may  be  represented  by 
error  boxes  as  shown  in  Figure  B2b.  An  infinite  number  of  straight  lines 
may  be  drawn  through  all  of  the  boxes.  The  lines  with  the  maximum  and 
minimum  slopes  are  shown  in  Figure  B2c  and  the  maximum  and  minimum 
intercepts  are  shown  in  Figure  B2d.  Thus,  the  maximum  and  minimum  for 
the  slope  and  intercept  are  determined  from  these  four  lines.  These 
"bands"  of  error  are  carried  along  through  the  rest  of  the  analysis  by 
using  them  in  the  next  calculation  or  plot. 

The  types  of  calculations  that  are  made  in  the  analysis  are  of  the 

form: 

Fave  * <1/4>  tl  F„  ‘“I 


1 
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F ■ G x H 

[B2] 

F « G ♦ H 

[B3] 

where  F,  G,  and  H are  variables. 

Equation  [Bl]  is  used  to  obtain  an  average  of  several  measurements 
of  the  same  thing.  In  the  analysis,  the  four  different  bias  moments 
were  used  to  obtain  (usually)  four  sets  of  data  for  each  velocity. 

Average  values  for  dC^/det  and  d^ /dec  are  obtained  and  used  in  Figure  3. 
In  that  figure,  the  uncertainty  in  the  average  value  is  assumed  to  be 
the  standard  deviation  of  the  four  (and  sometimes  fewer)  measurements. 

For  the  calculations  represented  by  Equations  [B2]  and  [B3] , 
respectively,  the  error  may  be  determined  from 

jdF/F | |dG/G | ♦ JdH/H J for  multiplicating  or  division 

jdF  | ^ | dG  | JdH  J for  addition  or  subtraction 

The  maximum  estimated  errors  determined  in  the  above  manner  for 
the  "indirect"  analysis  are  shown  in  Figures  B3  through  B6. 


From  Figures  B3  through  B6,  it  is  seen  that  the  estimated  maximum 
errors  are  considerable.  One  has  to  remember  that  these  are  maximum 
errors  possible.  The  major  cause  of  these  errors  is  attributed  to  the 
somewhat  inaccurate  measurement  of  absolute  flap  angle  (^).  By  tracing 
through  the  analysis  technique  shown  schematically  in  Figure  Bl,  it  is 
seen  that  this  variable  is  used  in  the  calculation  of  all  of  the  coef- 
ficients, centers  of  pressure,  and  effectivenesses.  Thus,  it  is  concluded 
that  in  order  to  use  this  analysis  technique  in  future  experiments,  it  is 
imperative  to  improve  the  measurement  accuracy  for  the  flap  angle  (while 
not  degrading  other  measurements) . 

As  can  be  seen  from  Figure  Bl,  the  foil  angle  (oO  and  the  velocity 
(U)  are  used  in  all  calculations  also.  Thus,  the  accuracy  of  these 
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MOMENT  CURVE  SLOPE  IN  RAD"  AND  MOMENT  EFFECTIVENESS 
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parameters  should  be  maintained  at  the  present  level,  or  perhaps  an 
improvement  in  the  foil  angle  (*)  would  be  in  order. 


The  accuracy  of  the  bias  moment  (M)  affects  the  moment  coefficients 
and  the  moment  effectiveness.  The  lift  (L)  affects  the  lift  coefficient 
and  the  lift  effectiveness.  The  centers  of  pressure  are  affected  by 
both  the  lift  and  moment  accuracies.  A possible  improvement  here  would 
be  to  improve  the  moment  accuracy  by  eliminating  its  variance  with  foil 
angle  (linkage  and  weight  position). 
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DTNSROC  ISSUES  THREE  TYPES  OF  REPORTS 

1.  DTNSRDC  REPORTS.  A FORMAL  SERIES,  CONTAIN  INFORMATION  OF  PERMANENT  TECH- 
NICAL  VALUE.  THEY  CARRY  A CONSECUTIVE  NUMERICAL  IDENTIFICATION  REGARDLESS  OF 
THEIR  CLASSIFICATION  OR  THE  ORIGINATING  DEPARTMENT. 

2.  DEPARTMENTAL  REPORTS,  A SEM1FORMAL  SERIES,  CONTAIN  INFORMATION  OF  A PRELIM- 
INARY, TEMPORARY,  OR  PROPRIETARY  NATURE  OR  OF  LIMITED  INTEREST  OR  SIGNIFICANCE. 
THEY  CARRY  A DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION. 

3.  TECHNICAL  MEMORANDA,  AN  INFORMAL  SERIES,  CONTAIN  TECHNICAL  DOCUMENTATION 
OF  LIMITED  USE  AND  INTEREST.  THEY  ARE  PRIMARILY  WORKING  PAPERS  INTENDED  FOR  IN 
TERNAL  USE.  THEY  CARRY  AN  IDENTIFYING  NUMBER  WHICH  INDICATES  THEIR  TYPE  AND  THE 
NUMERICAL  CODE  OF  THE  ORIGINATING  DEPARTMENT.  ANY  DISTRIBUTION  OUTSIDE  DTNSRDC 
MUST  BE  APPROVED  BY  THE  HEAD  OF  THE  ORIGINATING  DEPARTMENT  ON  A CASE-BY-CASE 
BASIS. 


